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ABSTRACT 

Taking into account the recent suggestion that a short gamma-ray burst 
(GRB) looks like the first 1 sec of a long GRB, we propose that the jet of a GRB 
consists of multiple sub-jets or sub-shells (i.e., an inhomogeneous jet model). The 
multiplicity of the sub-jets along a line of sight n s is an important parameter. If 
n s is large (^> 1) the event looks like a long GRB, while if n s is small (~ 1) the 
event looks like a short GRB. If our line of sight is off-axis to any sub-jets, the 
event looks like an X-ray flash or an X-ray rich GRB. The log-normal distribution 
of durations of short and long GRBs are also suggested in the same model. 

Subject headings: gamma rays: bursts — gamma rays: theory 

1. Introduction 

For the long gamma-ray bursts (GRBs), the cosmological distance, the collimated jet, 
the massive star progenitor, and the association with the supernova are almost established 
or strongly suggested (e.g., Meszaros 2002; Zhang &; Meszaros 2003). However, for short 
GRBs, little is known since no afterglow has been observed. The origin of the X-ray flashes 
(XRFs) also remains unclear although many models have been proposed (see Yamazaki, 
Ioka, & Nakamura 2004, and references therein). The observed event rate of short GRBs 
is about a third of the long GRBs while the observed event rate of XRFs is also about a 
third (Heise et al. 2001; Kippen et al. 2002; Lamb et al. 2003). Although there may be a 
possible bias effect to these statistics, in an astrophysical sense, these numbers are the same 
or comparable. If these three phenomena arise from essentially different origins, the similar 
number of events is just by chance. While if these three phenomena are related like Seyfert 1 
and 2 galaxies, the similar number of events is natural and the ratio of the event rate tells 
us something about the geometry of the central engine (Awaki et al. 1991; Antonucci 1993; 
Urry & Padovani 1995). In this Letter, we propose a unified model in which the central 
engine of short GRBs, long GRBs and XRFs is the same and the apparent differences come 
essentially from different viewing angles. 
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2. Unified Model 

It is suggested that short GRBs are similar to the first 1 sec of long GRBs (Ghirlanda, 
Ghisellini, & Celotti 2003). Although short GRBs are harder than long GRBs (Kouveliotou 
et al. 1993), this difference is mainly due to the difference in the low-energy spectral slope 
while the peak energy is similar (Ghirlanda, Ghisellini, & Celotti 2003). Other properties, 
such as (V/V max ), the angular distribution, the energy dependence of duration and the hard- 
to-soft spectral evolution of short GRBs, are also similar to those of long GRBs (Lamb et 
al. 2002). If short GRBs also obey the peak energy- luminosity relation found for the long 
GRBs (Yonetoku et al. 2003), it is suggested that short and long GRBs have a similar redshift 
distribution 1 (Ghirlanda, Ghisellini, & Celotti 2003). 

These similarities suggest that the difference between short and long GRBs is just the 
number of pulses, and each pulse is essentially the same (Ramirez- Ruiz & Fenimore 2000). 
As shown in Fig. 1, using 4Br catalogue of BATSE (Paciesas et al. 1999), the fluence is 
roughly in proportion to the duration in the range of 0.01 to 1000 sec (see also Balazs et 
al. 2003). Thus, we may consider that each pulse is produced by essentially the same unit 
or the sub-jet, and the GRB jet consists of many sub-jets. If many sub-jets point to our 
line of sight, the event looks like the long GRB while if a single sub-jet points to us, the 
event looks like a short GRB. Since we can observe only the angular size of ~ 7 -1 within the 
GRB jet with the Lorentz factor 7, different observers will see different number of sub-jets 
depending on the distribution of sub-jets within the GRB jet. Since the angular size of a 
causally connected region is also 7" 1 < 0.01, the opening half- angle of a sub-jet can be much 
smaller than that of the whole GRB jet (~ 0.1), say ~ 0.02. 

XRFs also appear to be related to GRBs. Softer and dimmer GRBs smoothly extend 
to the XRFs (Heise et al. 2001; Kippen et al. 2002; Lamb et al. 2003; Watoson et al. 2004), 
while the peak energy-isotropic luminosity/energy relations hold for GRBs as well as XRFs 
(Sakamoto et al. 2003; Yonetoku et al. 2003; Amati et al. 2002). The total energy including 
the radio afterglow of XRF 020903, which has a measured redshift, might be similar to that of 
GRBs (Soderberg et al. 2003). Other properties, such as the duration, the temporal structure 
and the Band spectrum of the XRFs are also similar to those of the GRBs, suggesting that 
XRFs are in fact soft and dim GRBs. In the sub-jet model, XRFs are naturally expected 
when our line of sight is off-axis to any sub-jets (Nakamura 2000; Ioka & Nakamura 2001; 
Yamazaki, Ioka, & Nakamura 2002, 2003b, 2004). 



1 Even if the afterglows of the short and long GRBs have a similar mechanism, the current limits are still 
consistent with the lack of afterglows for short GRBs (Hurley et al. 2002; Lamb et al. 2002; Klotz, Boer, & 
Atteia 2003). 
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The origin of sub-jets is not yet clear. In this Letter, we do not discuss the origin of the 
sub-jets, but argue the implications of the sub-jet model. 

3. An Example of Numerical Simulation of Our Unified Model 

We first show a numerical simulation to demonstrate how the event looks so different 
depending on the viewing angle in our unified model. Let us consider N tot = 350 sub-jets, 
for simplicity, confined in the whole GRB jet whose axis is the same as a d = axis. For 
each sub-jet the emission model is the same as in Yamazaki, Ioka, & Nakamura (2003b). 
Let the opening half-angle of the j-th sub-jet (j — 1, • • • , N tot ) be A#^ b , while the opening 
half-angle of the whole jet be A# tot . The direction of the observer and the axis of the 
j-th sub-jet are specified by (i? bs , y?obs) and <p U) ), respectively. We assume the j-th 
sub-jet departs at time t^ p from the central engine and emits at radius r = r (j> and time 
t = t U) = + r U) /(3 U) c, where t and r are measured in the central engine frame and we set 
^dep^ = 0- For simplicity, all sub-jets are assumed to have the same intrinsic properties, that 
is A6™ h = 0.02 rad, Y J) = 100 , r 0) = 10 14 cm, a$ = -1, (3$ = -2.5, ^hv'f = 500 keV and 
the amplitude A U) = const, for all j. The departure time of each sub-jet, t ( £. p is randomly 
distributed between t — and t = tdur, where tdur is the active time of the central engine 
measured in its own frame and set to t^ur — 30 sec. The opening half-angle of the whole 
jet is set to A6 to t = 0.2 rad as a typical value. We consider the case in which the angular 
distribution of sub-jets is given by P($ (j) , y? (j) ) d$ {i) d(p U) oc exp[-($ (3) /$c) 2 / 2 ] d$ {3) d(p U) for 

< A0 tot — A# sub , where we adopt $ c = 0.1 rad (Zhang et al. 2003). In this case, sub-jets 
are concentrated on the i? = axis (i.e., the multiplicity in the center n s ~ 10). For our 
adopted parameters, sub-jets are sparsely distributed in the range -# c < i? < A9 to t, however, 
the whole jet would be entirely filled if the sub-jets were uniformly distributed (i.e., the 
mean multiplicity n s ~ 3). Therefore, isolated sub-jets exist near the edge of the whole 
jet with the multiplicity n s <C 1 and there exists a viewing angle where no sub-jets are 
launched. Figures 2, 3 and 4 show the angular distributions of sub-jets and the directions 
of four selected lines of sight, the observed time- integrated spectra, and the observed light 
curves in the X-ray and 7-ray bands, respectively. Note here in Figure 2, "A" represents the 
center of the whole jet and is hidden by the lines of sub-jets. 

Long GRB: When we observe the source from the d = axis (case "A"), we see spiky 
temporal structures (Fig. 3- A) and E p ~ 300 keV which are typical for the long GRBs. We 
may identify case "A" as long GRBs. 

XRF and X-ray rich GRB: When the line of sight is away from any sub-jets (cases "Bi" 
and "B 2 "), soft and dim prompt emission, i.e. XRFs or X-ray rich GRBs are observed with 
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E p — 10 ~ 20 keV and ~ 4 orders of magnitude smaller fluence than that of case "A" (Fig. 
2). The burst duration is comparable to that in case "A". These are quite similar to the 
characteristics of XRFs. We may identify the cases "B^' and "B 2 " as XRFs or X-ray rich 
GRBs. 

Short GRB: If the line of sight is inside an isolated sub-jet (case "C"), its observed pulse 
duration is ~ 50 times smaller than case "A" (Fig. 3-C). Contributions to the observed 
light curve from the other sub-jets are negligible so that the fluence is about a hundredth 
of the case 'A" . These are quite similar to the characteristics of short GRBs. However the 
hardness ratio (= .5(100 — 300 keV)/S'(50 — 100 keV)) is about 3 which is smaller than the 
mean hardness of short GRBs (~ 6). Ghirlanda, Ghisellini, & Celotti (2003) suggested that 
the hardness of short GRBs is due to the large low-energy photon index a B ~ —0.58 so that 
if the central engine launches a B ~ —0.58 sub-jets to the periphery of the core where n s is 
small, we may identify the case "C" as the short-hard GRBs. In other words, the hardness 
of 3 comes from as = —1 in our simulation so that if as ~ —0.58, the hardness will be 6 
or so. We suggest here that not only the isotropic energy but also the photon index may 
depend on Another possibility is that if short GRBs are the first 1 sec of the activity of 
the central engine, the spectrum in the early time might be ~ —0.58 for both the sub-jets 
in the core and the envelope. This is consistent with a high KS test probability for E p and 
as (Ghirlanda, Ghisellini, & Celotti 2003). These possibilities may have something to do 
with the origin of as ~ — 1 for the long GRBs. 

X-ray pre/post- cursor. It is quite interesting that in Figure 4, we see the X-ray precursor 
at T bs ~ 60 sec in "B 2 " and the postcursor at T Q b s ~ 65-75 sec in "Bi". These can be 
understood by the model proposed by Nakamura (2000). 

4. Log-normal Distributions 

The total duration of the long and short GRBs are consistent with the log-normal dis- 
tributions (McBreen et al. 1994). In our sub-jet model, these distributions may be naturally 
expected as a result of the central limit theorem. Suppose a certain quantity q is expressed 
by a product of random variables q = x\Xi ■ ■ ■ x n . Then, log q = logxi + log x<i + • • • + logx n . 
When the individual distributions of log Xi satisfy certain weak conditions, the distribution 
of log q obeys the normal distribution in the limit of n — > oo by the central limit theorem. 
However in some cases the log-normal distributions can be achieved only by a few variables 
(Ioka & Nakamura 2002). Thus we might say, 'Astrophysically, not n — > oo but n = 3 
gives the log-normal distribution in practice!" . This argument may apply to the log-normal 
distributions of the peak energy, the pulse fluence and the pulse duration of GRBs (Ioka & 
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Nakamura 2002). 

In our sub-jet model, the short GRBs are due to a single sub-jet. The pulse duration 
of a single sub-jet is mainly determined by the angular spreading time scale and is given 
by the product of four variables in the internal shock model (Ioka & Nakamura 2002) as 
AT s hort ~ (1 + z) (L/c) (7 S /7 m ) 2 where L, 7 S , 7 m are the separation of two shells, the Lorentz 
factor of the slow and merged shell in the internal shock model, respectively. Therefore the 
log-normal distribution of the duration of the short GRBs may be a natural result of the 
central limit theorem. 

In our unified model, the duration of long GRBs is determined by the interval between 
pulses At = L/c times the multiplicity of the sub-jets n s . For a GRB at redshift z, the 
observed duration is given by the product of three random variables, ATi ong ~ (l+z)(L/c)n s . 
Therefore the log-normal distribution of the duration of long GRBs may be realized. The 
ratio of the duration of long GRBs to short GRBs is given by n s {j m /^ s ) 2 ~ 10 2 . Since in 
the internal shock model the relative Lorentz factor is not large, this equation suggests that 
n s = 10 ~ 30 which is compatible with the observed number of spikes of the long GRBs. 

5. Discussions 

Let A# su b, $ c and n s be the typical opening half-angle of the sub-jet, the core size of 
the whole jet and the mean multiplicity in the core. Then the total number of the sub-jets 
(N tot ) is estimated as N tot = n s (d c / A9 SVl \/) 2 ~ 10 3 so that the total energy of each sub-jet is 
~ 10 48 erg. In our model, the event rate of long GRBs is in proportion to $ 2 . Let M be the 
number of sub-jets in the envelope of the core with a small multiplicity n s < 1. Then the 
event rate of short GRBs is in proportion to MA9 2 uh so that M ~ 10 is enough to explain 
the event rate of short GRBs. 

Of course, the above numerical values are typical ones and should have a dispersion 
(Lloyd- Ronning, Dai, & Zhang 2003). Our core-envelope sub-jet model can have a similar 
structure to the two component jet model (Berger et al. 2003; Huang et al. 2003; Zhang, 
Woosley, & Heger 2003; Ramirez-Ruiz, Celotti, & Rees 2002) by varying such as n s and 
M. However the distribution of sub-jets could also have other possibilities, e.g., a hollow- 
cone distribution like a pulsar, a power law distribution, a Gaussian distribution (Zhang & 
Meszaros 2002; Rossi, Lazzati, & Rees 2002; Zhang et al. 2003) and so on. 

Some observers could see a cold spot with small n s in the core to have a small geomet- 
rically corrected energy even if the total energy of the GRBs is the same. Thus our model 
may be compatible with the recent claim that the total kinetic energy has smaller dispersion 
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than the geometrically corrected 7-ray energy (Berger et al. 2003; Bloom, Frail, & Kulkarni 
2003). The X-ray pre-/post-cursor is also expected if off-axis sub-jets are ejected earlier (for 
precursor) or later (for postcursor) than the main sub-jets (Nakamura 2000). The viewing 
angle of the sub-jets may also cause the luminosity-lag/variability/width relations of the 
GRBs including GRB 980425 (Yamazaki, Yonetoku, & Nakamura 2003; Ioka & Nakamura 
2001). This multiple sub-jet model is an extreme case of the inhomogeneous or patchy shell 
model (Kumar & Piran 2000; Nakamura 2000). The afterglow variabilities, such as in GRB 
021004, may arise from the angular energy fluctuations within the GRB jet (Nakar & Piran 
2003; Piran, Nakar, & Granot 2003), which might correspond to the inhmogeneous n s . 

Since the core may be regarded as a uniform jet, our model for the XRFs is analogous to 
the off-axis uniform jet model (Yamazaki, Ioka, & Nakamura 2002, 2003b, 2004). However 
the afterglow could have a different behavior between the core-envelope sub-jet model and 
the uniform jet model. In the uniform jet model, the afterglows of XRFs should resemble the 
orphan afterglows that initially have a rising light curve (e.g., Yamazaki, Ioka, & Nakamura 
2003a; Granot et al. 2002). An orphan afterglow may be actually observed in XRF 030723 
(Huang et al. 2003), but the light curve may peak too early (Zhang et al. 2003). The optical 
afterglow of XRF 020903 is not observed initially (< 0.9 days) but may not be consistent 
with the orphan afterglow (Soderberg et al. 2003). These problems could be overcome by 
introducing a Gaussian tail with a high Lorentz factor around the uniform jet (Zhang et 
al. 2003) since the energy redistribution effects may bring the rising light curve to earlier 
times (Zhang et al. 2003; Kumar & Granot 2003). The afterglow of a short GRB is difficult 
to predict since it could resemble both the orphan and normal afterglow depending on the 
sub-jet configuration within the envelope. 

Since all bursts have the same progenitor, our model suggests that the short GRBs and 
the XRFs are also associated with supernovae. The radio calorimetry will also give a similar 
energy to long GRBs because of the same reason. Our unified model will be refuted if the 
locations of short GRBs are mainly in the halo of the galaxy, as in the coalescing binary 
neutron star model (Bloom, Kulkarni, & Djorgovski 2002). 

Interestingly our model also predicts off-axis short GRBs or short XRFs. However these 
bursts will be difficult to detect since short XRFs, which have a multiplicity of n s ~ 1, 
will be ~ 30 times dimmer than XRFs with n s ~ 30. Note that the short XRFs will be 
longer than the short GRBs since the pulse duration grows as the viewing angle increases 
(Ioka & Nakamura 2001; Yamazaki, Ioka, & Nakamura 2002). The event rate of short XRFs 
will depend on the configuration of the sub-jets in the envelope. Further observations are 
necessary to determine the envelope structure. 
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Fig. 1. — The fluence S(50 — 300 keV) as a function of Tgo duration for BATSE bursts from 
4Br catalog. Courtesy of Drs. S. Michikoshi and T. Suyama. 
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Fig. 2. — The angular distribution of N tot = 350 sub-jets confined in the whole GRB jet 
in our simulation. The whole jet has the opening half-angle of A6 tot = 0.2 rad. The sub- 
jets have the same intrinsic luminosity, opening half-angles A# su b = 0.02 rad and other 
properties; 7 = 100, r = 10 14 cm, «g = — 1 , /3 b = —2.5 , hrfiS = 500 keV. The axes and the 
angular size of sub-jets are represented by crosses and the dotted circles, respectively. "A" 
represents the center of the whole jet and is hidden by the lines of sub-jets. 



- 12 - 




1 10 100 1000 
log [ hv (keV) ] 

Fig. 3. — Time-integrated energy spectrum of the emission from the multiple sub-jets for 
the observers denoted by "A", "B]", "B 2 ", and "C" in Figure 2. The source are located at 
z = l. 
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Fig. 4. — The observed X-ray and 7-ray light curves from the multiple sub-jets, corresponding 
the cases "A" (the upper left), 'TV' (the upper right), "B 2 "(the lower left) and "C"(the lower 
right) in Figure 2. The sources are located at z — 1. 



